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vn : 1. INTRODUCTION 



It has become common practice in nuclear physics to make a distinction between al- 
gebraic and geometric models. This is convenient for some purposes but, in general, 



(^1 it is misleading because almost all models in quantum mechanics have expressions in 

cn ■ algebraic terms. In fact, many-particle quantum mechanics is fundamentally a unitary 

representation of the Lie algebra of one-body operators and its observables are polyno- 
mials in the elements of this Lie algebra. 
K^ \ In proceeding towards a microscopic understanding of some nuclear phenomena, it is 

H ■ profitable to follow a sequence of steps along the following lines: 

(i) Observe the phenomena in many situations until its pattern becomes evident. 

(ii) Develop a phenomenological model to explain the observations and suggest new 
observations to ascertain the consistency of the model. 

(iii) Repeat steps (i) and (ii) to refine the model, assess its reliability and domain of 
validity, and identify its algebraic structure. 

(iv) If a model provides a successful understanding of the phenomena, the next step is 
to understand the model in terms of interacting nucleons. This is achieved if one is 
able to express the observables of the model in terms of nucleon coordinates and 
momenta. In the process, it is likely that the model will only be realisable in some 
limit or as an approximation to a more complex microscopic theory. It will then be 
of fundamental importance to determine if the limitations are consistent with the 
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FIGURE 1. The low-lying states of '°'*Er arranged into rotational bands 



experimental observations. As the following examples will illustrate, this exercise 
can also lead to improved phenomenological models with clearer microscopic 
foundations. 
(v) Having found a microscopically successful model, whose representations are capa- 
ble of explaining a range of physical observations, the next goal is to make use of 
its algebraic structure to identify an associated shell-model coupling scheme that 
will be appropriate for a shell-model description of the phenomena with realistic 
nucleon-nucleon interactions. 

The remarkable fact is that by such sequences of steps, one has obtained shell model 
coupling schemes appropriate for some of the dominant characteristics of nuclei. For ex- 
ample, the basic harmonic-oscillator shell-model, proposed to explain the strong bind- 
ing energy of so-called magic nuclei in terms of closed-shell structures, has provided 
the independent-particle basis for the standard shell model given by nucleons in a mean- 
field with spin-orbit interactions. Pairing models of singly closed-shell nuclei lead to the 
standard jj-coupling scheme of Flowers [1]; and nuclear rotational models lead to the 
Elliott SU(3) coupling scheme [2] for light nuclei. 

In the following we start from the observation of rotational bands organised as shown, 
for example, in Fig. 1 and the interpretations of such bands in the Bohr model and pro- 
ceed to their implications for a shell-model theory of strongly deformed nuclei. How- 
ever, before embarking on this process, we first recall some of the primary symmetry 
concepts that will be used. 



2. SYMMETRY GROUPS, DYNAMICAL GROUPS AND 
SPECTRUM GENERATING ALGEBRAS 

Definition: A symmetry group of a system is a group of transformations of the system 
that leave its Hamiltonian invariant. 

There are several variations in the literature on the definition of a dynamical group 
[3]. The following definition is most useful for present purposes. 

Definition: A dynamical group for a Hamiltonian H is a Lie group of unitary transfor- 
mations of the Hilbert space H of a system with Hamiltonian H such that the subspaces 
of EI that are invariant under the transformations of the dynamical group are spanned by 
eigenstates oiH. 

Consider, for example, a particle in ordinary three space moving in a central force 
potential V{r) such that the Hamiltonian is S0(3) invariant. Basis wave functions for 
the particle's Hilbert space are given by products of radial wave functions and spherical 
harmonics 

Wnim{r, 0, (p) = Rni{r)Yi„{9(p), (1) 

where / and m are angular-momentum quantum numbers and n is a radial quantum 
number. The Hamiltonian matrix in this basis is given by 

HL,n'm' = {¥nlm\H\\ifn>Un') (2) 

which, because S0(3) is a symmetry group, has Z as a good quantum number. In addition, 
because H is S0(3) invariant, the Hamiltonian matrix elements are diagonal in m and 
independent of its value, i.e., 

and m is also a good quantum number. Thus, to obtain the energy spectrum for the 
particle, it remains only to diagonalise each //' with respect to the radial quantum 
number n to obtain its eigenfunctions 

^alm = / .^anWnlm- (4) 

n 

This process is greatly simplified by the fact that the central force problem has a 
dynamical group SU(1,1) for which the set of states {v^„/m;n = 0, 1, 2, . . .} for each pair 
of Im values carries an irreducible representation. 

Definition: A spectrum generating algebra (SGA) for a Hamiltonian is a Lie algebra 
such that the Hamiltonian can be expressed (most usefully as polynomials) in terms of 
its elements. 

An SGA Q is most useful because its representation theory enables one to calculate 
the matrix elements of any X E Q and hence the matrix elements of the Hamiltonian. A 



Hamiltonian may have many SGA's. A particularly useful choice is one for which the 
important observables of the system have simple expressions in terms of its Lie algebra 
Q, e.g., as linear or quadratic polynomials of Lie algebra elements. 

An SGA for the above central-force problem is the Lie algebra su(l,l) whose in- 
finitesimal generators include the square of the radial coordinate, r^, of the particle and 
the square of it momentum, p^. 

A system may have many symmetry groups and many dynamical groups. For exam- 
ple, the group generated by an SGA for a Hamiltonian is a dynamical group for that 
Hamiltonian. 



3. THE BOHR COLLECTIVE MODEL 

The Bohr model is a liquid-drop model with quadrupole shape coordinates defined for 
its surface radius by 

R{e,(p)=Ro\l+Y,cCyYiv{d,(p) + ..- 

V 

It is quantised by the introduction of quadrupole shape observables {dv} and canonical 
momenta {ft^} that satisfy the Heisenberg commutation relations 

and act as operators on a Hilbert space HI of square-integrable wave functions of the 
shape variables, according to the equations 

dvWa) = avW«), ft^ = -ih- — v/(a), v = 0, ±1,±2. 

oav 

A Hamiltonian for the model is then an S0(3)-invariant operator of the form 

B 

where V{a) is a rotationally-invariant potential energy (a function of the coordinates) 
and 5 is a mass parameter. 

The model has geometric and algebraic structures that make it easy to use, but which 
have only recently been exploited in what is known as the Algebraic Collective Model 
(ACM) [4, 5, 6, 7]. The coordinate space of this model is the real five-dimensional 
Euclidean space M^. This space has much in common with the more familiar three- 
dimensional space R-^ and similarly, for a rotationally-invariant Hamiltonian, it is most 
naturally assigned spherical polar coordinates. Apart from the origin (which is a point 
of measure zero), the Euclidean space R^ is a tensor product of a radial line and a two- 
sphere of unit radius as illustrated in Fig. 2(a). Similarly, to within a point of measure 
zero, the Euclidean space R^ is a tensor product of a radial line and a four-sphere of 
unit radius. Thus, whereas R^ is assigned spherical polar coordinates (r, 0,^), where 
r is a radial coordinate and (0,(p) are coordinates for a point on the two-sphere, the 
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FIGURE 2. Spherical polar coordinates (a) for a three-dimensional and (b) for a five-dimensional 
Euclidean space. 



space M^ is assigned spherical polar coordinates (/3,7, Q), where now /3 is a radial 
coordinate, defined by /3^ = a ■ a, which together with an angle coordinate 7 determine 
the triaxial shape of a rotor; Q. denotes the orientation angles of this rotor. Together, the 
angle coordinates (7, Q.) complete a system of spherical coordinates for the four-sphere. 
Thus, whereas the Hilbert space for IS? is the tensor product 



^2(^3) ^^2 



X if 2 (52) 



(5) 



and is spanned by functions {v^,j/ra(r, 0, (p) = Rin{r)Yi^{Q , (p)}, where {7?/,,(r)} is abasis 
for a unitary irrep (irreducible representation) of the group SU(1,1) and {y/,„(0, 9)}is a 
basis of spherical harmonics for the two-sphere, the Hilbert space for M^ is the tensor 
product 

if2(^5^) _ ^2(^^^) ^ ^2(5^) (6) 

and is spanned by functions {Mvn{P)^vaLM{Yi^)} , where {=^v«(j8)} is again a basis 
for a unitary irrep of the group SU(1,1) and {^vaLMij,^)} is abasis of S0(5) spherical 
harmonics for the four-sphere. 

Matrix elements of operators of interest on the Hilbert space ^^(R-') can be calcu- 
lated algebraically. This is because explicit expressions are known [8] for the matrix 
elements of r, 1/r, r^, l/r^, V^, for S0(3) spherical harmonics, for Clebsch-Gordan co- 
efficients, and for the S0(3)-reduced matrix elements 



(YlMWYh) 



(2Z2 + l)(2Zi + i; 



4n 



(/1O/2OI/3O). 



(7) 



Similar expressions are now known for the Hilbert space Jf^{M.^) of the Bohr model. 
By simple algebraic methods [8], we obtain explicit expressions for the matrix elements 
of /3, l//3,/32, 1//32 and V^. From the basic S0(5) spherical harmonics 



^1122(7,^) - cos 7^22 + V2sinr(^l2(^) + ^\2(^)) 



(8) 



^2122(7,^) - cos2r^22 - V 5 sin2r (^|2(^) + ^^2 2(a)) , (9) 

^3ioo(r,^)-cos3r, (10) 

^3133(7,^)- A/isin3r(^|3(a)-^2^,3(a)), (11) 



we can generate the complete set of S0(5) spherical haraionics and, from them, derive 
all required S0(5) Clebsch-Gordan coefficients [9, 5]. We also obtain [5] the S0(5)- 
reduced matrix elements 



m\\\^\W) = 1 (f + ^)' (2vi+3)(2v2 + 3) 



, (a + 4)(a-2vi + l)!(a-2v2 + l)!(g-2v3 + l)! ,.^, 
""'^ {^^\ ' ^^^^ 

where o = v\ +V2 + V3. Thus, Bohr model calculations are carried out quickly and easily. 
Particularly informative are Bohr Hamiltonians are of the form 

^(5,A,z,'c) = -^ + y;t(/3,7) (13) 

with 

yx{^.y) = ^5[(1 -2A)/32 + A/34] -;t/3cos3r+Kcos2 3r. (14) 

For example. Fig. 3(a) shows the low energy-level spectrum and E2 transition rates for 
such a Hamiltonian with k = 0. The physical content of this spectrum is exposed by 
comparing it with results for the simpler Hamiltonian with potential 

y(i8,r)^V(/3o,ro) + ^5«2(/3_/3o)2 + l5/32£oJ(r-7o)', (15) 

where j8o = and 7o = are the values of j8 and y at the potential minimum. With the 
additional assumption that the rotational motions of the model are slow relative to its /3 
and /vibrational modes there is an adiabatic decoupling of the rotational and vibrational 
degrees of freedom and the spectrum that emerges becomes as shown in Fig. 3(b). Thus, 
the relatively small differences between Figs. 3(a) and 3(b) can be attributed to S0(5) 
centrifugal coupling effects between the rotational and vibrational degrees of freedom 
that are included in Fig. 3(a) but not in 3(b). 

The Bohr model and its many generalisations by Bohr, Mottelson, and numerous 
colleagues [10] has been enormously successful in the interpretation of a huge body 
of experimental data on collective motions. A particularly significant generalisation was 
the extension of the Bohr model to describe the dynamics of a nucleon coupled to an 
even-mass rotational nucleus [11]. This and other generalisations reveal how nucleons 
in a deformed mean field can form a rotor with many intrinsic degrees of freedom 
that can rotate adiabatically without undue disturbance of its intrinsic structure. Such 
generalisations are known as unified models. 
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FIGURE 3. (a) The low-energy spectrum of the HamOtonian H{B,X,x, K = 0) of Equation (14) for 
B = 20, /I = 1.5, and x = ±2.0. Reduced E2 transition rates are shown in units for which B(E2;2i — > 
Oi) = 100. Energy levels are given in units such that the lowest L~2 state has energy £21 = 6. (b) The 
corresponding results for the adiabatic Bohr model in the axially symmetric limit showing a ground-state 
band, a one-phonon beta vibrational band, and one- and two-phonon gamma vibrational bands. (Figure 
from [6]) 



4. PROGRESSION TO A MICROSCOPIC COLLECTIVE MODEL 

In view of the successes of the Bohr and unified models, a great deal of effort has been 
expended in seeking to express collective dynamics in terms of microscopic nucleon 
coordinates. In retrospect, this is simple if one starts by replacing the surface shape 
coordinates {ttv} of the Bohr model by microscopic Cartesian monopole/quadrupole 
moments 



{Jij 7 , ^ni^nit hj '^ ■> -^ ■> •-' ■> ^ ^^5 



(16) 



n=\ 



Time derivatives and corresponding momentum observables are then given by 
dQ,j 



Q, 



dt 



7 , [.^ni^nj ~r ^ni-^nj ) 7 ^ij "^ idij 7 , yPni^nj ~r ^niPiij J 7 K'^ ' ) 



where M is the nucleon mass. These shape and momentum observables are now micro- 
scopic and close on a Lie algebra with commutation relations. 



[Qij,Pki] = in [5iiQjk + 5ikQji + 5jiQ,k + 5jkQ,i 



(18) 



This Lie algebra differs from that of the Bohr model. However, the subset of 
quadrupole moments defined in a spherical tensor basis by 

^o = ;7|e(2eii-<222-<233), 

4±i = Vle(ei2±iei2), (19) 

^±2 = ^ £(222 -233 ±21223), 

have commutation relations with the corresponding quadrupole momenta that contract to 
those of the Bohr model when restricted to states of the Hilbert space of large monopole 
and relatively small quadrupole moments. For example, with the monopole moment 
defined by 

50to = <2ll + <222 + <233, (20) 

and £2 = l/(Sato), it follows that 

[q^,p'']=m^+0{£'^)^ih5^, ase^O. (21) 

We refer the model with this SGA as the microscopic Bohr Model. 

In spite of this result, this model fails to describe the low-energy rotational bands for 
which the Bohr model was introduced. This is because its momentum operators are the 
infinitesimal generators of irrotational-flow rotations for which the moments of inertia 
are much smaller than those needed to describe observed rotational bands. However, this 
limitation is overcome by enlarging the Lie algebra, spanned by {Qij^Pij}, to include the 
angular momentum operators 

Lfc = hLij = Y, {XniPnj - XnjPni) ■ (22) 

n 

The Lie algebra obtained is then isomorphic to a known CM(3) Lie algebra introduced 
by Weaver et al. [12]. 

The CM(3) model is richer than the Bohr model. Its six independent {Qij} observ- 
ables characterize the size and quadrupole shape of a nucleus and, among its 9 momenta 
{Lij.Pij}, the angular momenta Lij are infinitesimal generators of rigid-body rotations, 
the three components of Pij with i ^ j are infinitesimal generators of irrotational-flow 
rotations, and the three components Pa are infinitesimal generators of diagonal shape 
scale transformations [13]. Combinations of these momentum operators also generate 
intrinsic current circulations as illustrated in the third row of Fig. 4. 

The unitary irreps of CM(3) were determined by Rosensteel [14] with the intention 
of developing a microscopic theory of collective motion. However, two obstacles were 
encountered: one was that, while the many-nucleon kinetic energy has a simple known 
expression, 7" = 2^ UniPfit^ its CM(3) collective model component (derived for example 
in [15]) proved to be excessively complicated; another was that the CM(3) representa- 
tions did not relate in a natural way to the nuclear shell model. The way around both 
obstacles was simply to augment the CM(3) algebra further so that it included the many- 
nucleon kinetic energy. The result was the much more powerful symplectic model. 
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FIGURE 4. Collective flows generated by CM(3) transformations 



5. THE SYMPLECTIC MODEL 

The symplectic model [16, 17] is an algebraic model whose dynamical group Sp(3,R) 
comprises all n-independent linear canonical transformations of the single-particle 
phase-space observables 

J J 

that preserve the commutation relations [x„i,Pmj\ = ih5m.n8i.j- Its spectrum generating 
algebra is spanned by the elements of the CM(3) Lie algebra 



SJ-ij 7 , ^ni^njt ^ij 7 , \ Pni^n ;' ~l~ ^niPn ;' ) ; ^ ij 7 , [-^niPn j ^njPniJt \^^) 



n=\ 



plus the bilinear momentum operators 

^'j = Z^PniPnj- 



(24) 



Adding the Ktj operators to the CM(3) algebra is invaluable for many reasons: 

(i) The representations of Sp(3, R) are easier to handle than those of CM(3). 
(ii) The many-nucleon kinetic energy and harmonic-oscillator Hamiltonians 
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(25) 



where M is the nucleon mass, are elements of the Sp(3, R) Lie algebra for all values 
of coi, 0)2, ft)3. As a result, the symplectic model is compatible with the harmonic 
oscillator shell model and with unified collective models such as the Nilsson model, 
(iii) Sp(3,R) representation theory makes it possible to diagonalize a collective model 
Hamiltonian 

^ = 2^1^" +^(2)' (26) 

i 

in a basis of shell-model states. 

(iv) The group, Sp(3, R), contains important subgroup chains 

Sp(3, R) D CM(3) D R0T(3) D S0(3), 
Sp(3, R) D U(3) D SU(3) D S0(3), 

where R0T(3) is a dynamical group for a rigid-rotor model and SU(3) is the 
dynamical group of Elliott's SU(3) model. The first chain reflects its content in 
terms of the dynamical flows of a quantum fluid. The latter chain defines an optimal 
shell-model coupling scheme for a shell-model description of collective states in 
nuclei. 

These properties are all you could wish for from a microscopic collective model. 



5.1. The Sp(3, R) algebra in a U(3) basis 

The U(3) C Sp(3, R) subalgebra consists of all elements of the Sp(3,M) Lie algebra 
that commute with the spherical harmonic-oscillator Hamiltonian 



MA) _ Y ^ 



l^^l+^-M(D\xl,+xl^+xl,) 



(27) 



where co is fixed to get the mean-square radius of a nucleus at its observed value. 

To identify this subalgebra, express the Xni and ^,„- observables in terms of harmonic 
oscillator raising and lowering operators 

-Wni + Cni), Pni = ih^={clf-Cni), (28) 



V2a"" ' "'^' "" V2' 



where a = ^/Wcojh is the harmonic oscillator unit of inverse length. These operators 
satisfy the commutation relations 

k"/, clj] = d,n,„dij, [h'^^o^ cl^ = hcocl-, [H^^Q, c,„] = -hcocni. (29) 

Expansion of the Sp(3,i?) operators in terms of them gives 

a^Q,j = Aj+h + Bij. (30) 



Pij = 2ifi{Aij-Bij), 
where 

^ij = L?i \Snfni +27' ~^0' ^ 2 V^O' + ^ji) ' 

satisfy the commutation relations 

ruW J. .1-0 \H^^^ L--1-0 



(31) 



(32) 



Thus, the U(3) subalgebra is spanned by the subset of operators {.^,j, L//} that commute 

with H^yIq, and the operators {Aij} and |5,j| are, respectively, ±2ho} raising and 
lowering operators for representations of the sp(3,]R)) algebra. The L = combination 

Aq = y |L,A,/ is the creation operator for a 2h(0 giant-monopole excitation and the 

five L = 2 combinations {A2v',y = 0, ±1,±2} are creation operators for 2hco giant- 
quadrupole excitations. Together, these six raising operators transform as the L = and 
2 components of an SU(3) (2 0) tensor (see Exercise ??). 



5.2. Irreducible representations (irreps) of Sp(3, R) in a U(3) basis 

The Sp(3, R) Lie alglebra is non-compact and, like the Lie algebra of the Bohr model, 
it has only infinite-dimensional unitary irreps. In the space of the nuclear shell model it 
has irreps with lowest weights, but none with highest weights. Thus, irreps of Sp(3,R) 
are constructed within the shell-model space by first defining their lowest- weight states 
which are also conveniently defined as U(3) highest-weight states. 

A U(3) highest- weight state is labelled by three quantum numbers |A^(A/i)) and 
satisfies the equations 

Cu\NiX^))=Ni\N{X^)), /= 1,2,3, (33) 

iCn+C22 + C33)\NiX^)) =N\N{X^)), (34) 

(Cii-C22)|A^(Aju))=A|^(Am)),, (35) 

{C22-C33)\N{X^)) = H\N{kn)), (36) 

Cij\N{?i^))=0, Vl</<j<3, (37) 

with 

A^ = A^i+A^2 + A^3, X=Ni-N2, ll=N2-N3. (38) 

An Sp(3,R) lowest-weight state \No{A{)iJ{))) is a U(3) highest-weight state that is also 
annihilated by the giant-resonance lowering operators, i.e., 

B,j\No{Ao^)) =0, VI < / < J < 3. (39) 



A U(3) irrep with highest- weight state \N{Xil)) is denoted by {N{Xil)} and an Sp(3, R) 
irrep with lowest-weight state |A''o(Ao/io)) is denoted by (A^ol-^oMo))- A U(3) irrep whose 
highest-weight state is also an Sp(3,R) lowest-weight state is said to be a lowest-grade 
U(3) irrep for that Sp(3,R) irrep. It follows that basis states for an Sp(3,R) irrep are 
obtained by augmenting a basis for a lowest-grade U(3) irrep by addition of an infinite 
set of giant-resonance excitations generated by the {Ajj} raising operators. 

There are four classes of U(3) irreps: {^(00)}, {N{XO)}, {A^(OA)} and {N{X[i)}. 
The {//(OO)} irreps are one-dimensional and contain a single L = state. The {N{OX)} 
and {A'^(AO)} irreps are contragredient to one another; they have common spectra of 
angular momentum states, given by 



L = A, A-2, A-4, ...,lorO, 



(40) 



and differ only in that their quadrupole matrix elements are of opposite sign. A generic 
U(3) irrep {N{Xil)} contains a single SU(3) irrep (A/i) and has a spectrum of S0(3) 
angular-momentum states given by the SU(3) — >■ SO (3) branching rule 



(Am) 



^ L = A, A-2, .. 


., lorO, wheni^ = 


h-> L = K,K+\, . 


.,K + X, wheni^^O 


for K = ji, 


/i-2, /i-4, ..., 1 orO. 



(41) 



At this point it is useful to recall that SU(3) irreps closely resemble truncated irreps 
of a rotor In fact, for values of A and/or /i large compared to the angular momenta of 
the states of interest, the properties of the SU(3) irreps (A 0) and (OA) approach those of 
prolate and oblate rigid rotors, respectively, whereas the properties of a generic SU(3) 
irrep (A/i) approach those of a triaxial rigid rotor. Thus, an SU(3) model was introduced 
by Elliott [2] as a first step towards a shell-model description of rotational states in light 
nuclei. Typical spectra for SU(3) model irreps are shown in Fig, 5. 
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FIGURE 5. Typical SU(3) model spectra for the lowest-energy states of three irreps of even values of 
Aq. The (0 0) irrep contains a single state. 

The energy spectrum of basis states for an Sp(3,R) irrep, with respect to the spheri- 
cal harmonic-oscillator Hamiltonian, Eq. (27), is now obtained as follows. The states of 



a lowest- grade U(3) irrep {{No{Xo iAq))} all have the common harmonic-oscillator en- 
ergy Noho). One-phonon giant-resonance states appear at an energy of {No + 2)h(X) and 
because the giant-resonance raising operators are components of an SU(3) (20) tensor, 
they generate the states of a generally reducible SU(3) representation given by the SU(3) 
tensor product (Ao/io) CS (20). The spectra of basis states are illustrated in Fig. 6. 
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FIGURE 6. Low basis states for Sp(3 , M) irreps ordered by their spherical harmonic oscillator energies: 
(a) for an {Nq{00)) irrep; (b) for a {No{XoO)) irrep; (c) for a generic {Nq{?lqI1o)) irrep. 

It is of interest to note that the spectrum of states for an (A'o(OO)) irrep is in 1-1 
correspondence with that of a 6-dimensional harmonic oscillator. In fact, they are a basis 
for an irrep of the microscopic Bohr model, spanned by the {Qij^Ptj} operators of Eq. 
(17), which comprises giant monopole and quadrupole harmonic vibrational states with 
irrotational-flow mass parameters. More generally, the states of an Sp(3,R) irrep are in 
1-1 correspondence with the states of such a Bohr model coupled to a non-trivial SU(3) 
irrep. This is illustrated pictorially in Fig. 7. Moreover, it has been shown [18, 19, 20] 
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FIGURE 7. Basis states for the Sp(3, R) model are in 1-1 corresondence with a coupled product of the 
Bohr model (realised in terms of giant monopole-quadrupole resonance states with irrotational flows) and 
an SU(3) model. 



that, in the limit of large quantum numbers, the Sp(3, R) model contracts to the coupled 
product of a Bohr vibrational model and a rotor model. This is the so-called coupled 
rotor-vibrator limit of the symplectic model. 

A simple prediction that follows, without a calculation, from the observation that the 
underlying microscopic description of monopole/quadrupole rotational and vibrational 
dynamics is provided by the symplectic model, is that this model contains no low-energy 
/3 or 7 vibrational bands [21]. It does exhibit excited K = 2 bands but these are more of 
a triaxial rotor than of a vibrational nature. Excited K — bands (for even-even) nuclei 
do occur in the model but are associated with an excited Sp(3, R) irrep. With mixing of 
Sp(3,R) irreps, such an excited band can exhibit properties associated with a /3 band. 
However, the occurrence of such bands is much more appropriately described in terms 
of band mixing and as a departure from the underlying symmetry rather than as a part 
of it. The study of such ^ = bands and the E2 transitions between them is therefore 
expected to provide valuable information on the mixing of Sp(3, R) irreps in even-even 
nuclei (see, e.g., [22]). 



5.3. A model Sp(3, R) calculation 

Figure 8 shows the energy levels and reduced E2 transition rates for a symplectic 
model calculation [23] for an Sp(3,R) irrep (826.5(78,0)) with the two-parameter 
Hamiltonian 

ff=<o'+z(e-e+2^). (42) 

The irrep (826.5(78, 0)) was chosen to be appropriate for ^^^Er on the basis of experi- 
mental data [24] (as discussed below). The two parameters £ and x of this Hamiltonian 
are naturally defined as follows: £ can be adjusted so that the potential energy compo- 
nent of the Hamiltonian (42) has a minimum at the deformation determined from ex- 
perimental E2 transition rates; the strength parameter x is then determined by the self- 
consistency condition that the density distribution for the wave function has the same 
quadrupole shape as that of the potential. For the results shown, the simpler procedure 
was followed of adjusting £ as mentioned above to accord with the observed deforma- 
tion, and adjusting x to give the observed energy of the 2i state. The close agreement of 
the results, especially for states of lower angular momentum is very satisfying because it 
means that the symplectic model is describing the physics of a rotational band of states 
with realistic moments of inertia that emerge naturally without any adjustment of the 
many-nucleon kinetic energy, e.g., by introducing an effective nucleon mass. 

The fact that the symplectic model agrees closely with a two-parameter rotor model is 
also satisfying because, while the rotational character of the states described is presumed 
in the rotor model, no such assumption is made in the symplectic model which has 
many more degrees of freedom. In fact, it is apparent that the experimental energy levels 
exhibit a degree of centrifugal stretching that is not admitted in the rigid-rotor model. 
Thus, the rotor model interpretation of the symplectic model results suggests that a better 
agreement with the higher energy levels could have been obtained in the symplectic 
model with a potential energy having a less rigidly-defined minimum. The agreement 
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FIGURE 8. Low-energy states of ^''''Er calculated in an Sp(3, R) irrep (826.5(78,0)) for the Hamilto- 
nian (42). For comparison the spectra calculated in the rotor and SU(3) models are also shown. The boxed 
numbers are B(E2)'s in W.u. 



between the two-parameter SU(3) and rigid-rotor models is not surprising, in view of the 
contraction, noted above, of the SU(3) model to the rotor model for large-dimensional 
representations. However, the relationship between the symplectic model and these two 
models is important for understanding the microscopic structure of rotational states and 
for understanding how the seeds of nuclear rotations can be found within the framework 
of a shell-model coupling scheme. 

6. THE SHELL-MODEL THEORY OF DEFORMED NUCLEI 



Because the many-nucleon Hilbert space is infinite dimensional, many-nucleon calcu- 
lations must necessarily be restricted to finite-dimensional subspaces. Thus, to obtain 
meaningful results, it is necessary to choose the finite-dimensional subspaces such that 
they have significant relationships with the physical states they purport to describe. 

The standard shell model is based on the presumption that the dominant forces on 
a nucleon in the nucleus are generated by a spherically symmetric mean-field with 
spin-orbit interactions. Thus, a hierarchy of standard shell model states is defined by 
a partially-ordered basis of independent-particle model states. In this model, many- 
nucleon subspaces are naturally defined by specification of occupied and unoccupied 
single-particle states and a complementary valence space of partially-occupied states. 
Such subspaces can be ordered by their summed independent-particle-model energies. 
Which subspace one chooses for a specific purpose obviously depends on one's objec- 
tives and computational resources. However, a common practice is to choose a subspace 



that is as large as can be handled and, apart from symmetry constraints (such as parity), 
is one for which all independent-particle states of energy below those of the valence 
space are occupied. 

In practice, the mean field of the shell model is commonly replaced by a spherical 
harmonic oscillator potential. This makes it possible to capitalise on the rich variety 
of symmetries associated with the harmonic oscillator, which greatly facilitate the exe- 
cution of complex calculations. It must be remembered though that the purpose of the 
independent-particle potential is to define basis states for a many-nucleon calculation 
and the better the basis the fewer basis states are needed to obtain results to a given level 
of accuracy. Except for the radial tails of the wave functions, which are important in re- 
action theory but less important for the calculation of bound states, a spherical harmonic 
oscillator does a reasonably good job of approximating the mean field for spherical nu- 
clei and the inadequacies of the bases that it provides are more than compensated by the 
many benefits resulting from the harmonic oscillator symmetries. 

A concern addressed in these lectures is that, while the standard spherical harmonic 
oscillator shell model is very successful for the description of spherical nuclei, it fails 
badly in the description of strongly-deformed nuclei. This is revealed by considering 
which Sp(3, R) irreps are appropriate for the description of rotational bands in nuclei 
and finding that they lie in shell-model subspaces that are often far above those that 
would normally be included in the valence spaces of the spherical harmonic-oscillator 
shell model. It is insightful to consider this concern from a mean-field perspective. 



6.1. Geometrical and mean-field perspectives 

Loosely speaking, mean-field theory can be regarded as quantum mechanics con- 
strained to an orbit of a dynamical group that is isomorphic to a classical phase space. 
Thus, mean-field theory provides an interface between classical and quantum mechanics 
and enables one to obtain valuable classical insights into the quantum world. For exam- 
ple, in quantum mechincs, the ground state of a Hamiltonian is contained in an irrep of 
a dynamical group for that Hamiltonian. In mean-field theory, a classical ground state is 
an equilibrium state given by the minimal-energy state on the classical phase space. This 
minimum energy is then an upper bound on the quantum-mechanical ground state en- 
ergy. Thus, it is meaningful to order Sp(3, R) irreps by their minimal energy mean-field 
states. 

Let I a) = lA^o('^oMo)) denote a lowest- weight state for an Sp(3,R) irrep and let t{g) 
denote a unitary transformation of the states of this irrep by an element g E Sp(3, R). 
The group orbit is then the manifold of coherent states 

^ = {\a{g)) = f{g)\a); g e Sp(3, R)}, (43) 

which is isomorphic to the classical phase space for the irrep. We now claim that the real 
submanifold of zero-momentum states of this phase space is the set 

M = {\o{g)) = t{g)\ay, g e GL+(3,R)}, (44) 



where GL+(3, R) is the connected subgroup of real 3x3 matrices with positive determi- 
nant. To see this, first observe that the phase space ^ has dimension dim(^) = 18; this 
follows because the Sp(3, R) Lie algebra is of dimension 21 but, as shown by Eqs. (34) - 
(36), three elements of this Lie algebra, namely C\\,C2i and Qs, leave the state |c7) in- 
variant. Next observe that the nine linearly-independent momentum operators Pij = Pji 



and L 



u 



Lji are the infinitesimal generators of the real subspace 



C 



and that 



they span a Lie algebra isomorphic to that of GL-|-(3, R). 

We now make use of the fact that a general linear matrix g G GL_|_(3,R) can be 
expressed as a matrix product 

g = Q.diQ!, (45) 

where ^ and Ql are S0(3) matrices and d = diag(Ji, ^2,^3) is a real diagonal matrix 
with positive entries. Thus, in the unitary representation T , we can express the mean- 
field states of ^ in the form 



\o['^dQ!))=R{Q.)U{d)R{Q!)\o), 



(46) 



where ^(Q) is a rotation/inversion operator and f/(d) is a unitary representation of the 
scale transformation {xni — ^ diXni, i= 1,2,3}. 

The sequence of transformations R{^)U {A)R{'D!) of the Sp(3,R) lowest-weight state 
\o) is represented pictorially in Fig. 9 for two Sp(3,R) irreps, (A'o(OO)) and a generic 
irrep (A^o('^Mo))- For the purpose of this representation, it is convenient to make use 
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FIGURE 9. A geometrical perspective of the nine dynamical degrees of freedom of the symplectic 
model which reduce to six for a spherical irrep {No{OQ)) wherein the three SU(3) rotational degrees of 
freedom are suppressed; see the text for an explanation of this figure. 

of the 1-1 correspondence between the basis states of the symplectic model and those 
of a coupled Bohr (giant-resonance) model and an SU(3) model as depicted in Fig. 
7. The states of the manifold ,^ can then be represented by a pair of ellipsoids: one of 
fixed shape but arbitrary orientation, representing the three rotational degrees of freedom 
of the SU(3) rotor; and one with a full range of ellipsoidal shapes and orientations, 
representing the six degrees of freedom of the Bohr model. The states of the SU(3) rotor 
are drawn in the figure with solid lines while those for the giant-resonance are drawn 



with dashed lines. Recall, however, that for the Sp(3, R) irrep {No{00)) the lowest-grade 
SU(3) irrep is the trivial identity representation. For this irrep only the giant-resonance 
degrees of freedom survive and the lowest-weight state (A^o(OO)) is the giant-resonance 
vacuum state which has a spherical density distribution as shown in row (i) column (a) 
of the figure. For a generic Sp(3,R) irrep the lowest- weight state |a) is again a giant- 
resonance vacuum state. But it is also the highest-weight state for a non-trivial SU(3) 
irrep with a non-spherical density distribution, as shown in row (ii) column (a). 

Because the action of the rotation operator R{Q.') on the lowest- weight state |(7) does 
not excite any giant-resonance states, its effect is simply to reorient the SU(3) ellipsoidal 
density distribution. Conversely, the transformation f/(d), shown in (c), is a diagonal 
scale transformation which excites giant resonance degrees of freedom and changes the 
spherical shape of the giant-resonance vacuum state to a more general ellipsoidal shape, 
with its principal axes aligned along the space-fixed axes, but leaves the SU(3) degrees 
of freedom alone. Finally the last rotation R{Q.), shown in (d) rotates both the SU(3) and 
giant-resonance ellipsoids. The resulting shape of the nucleus, as represented by Fig. 
9(d), is now a composition of the two ellipsoids shown in which that corresponding to 
the SU(3) lowest grade ellipsoid is scaled in a manner defined by the GR ellipsoid. 



6.2. The minimal-energy mean-field state for an Sp(3, R) irrep 

To derive the minimal-energy mean-field state for an Sp(3, R) irrep, one would ideally 
make a variational calculation, analogous to a Hartree-Fock calculation, on the phase 
space of Sp(3,R) coherent states defined by Eq. (43). However, it is possible to assess 
what the result would be on the basis of a simple self-consistency argument. First 
observe that, for a time-reversal-invariant Hamiltonian, the minimal-energy state should 
be one with vanishing expectation values of its momentum operators. Thus, it should be 
the state in the set {|a(^)), ^ e GL+(3,R)} for which 

E{ada') = {o{adQ.')\H\o{adQ.')) (47) 

is a minimum for the given Hamiltonian. 

The next observation is that, for a rotationally-invariant Hamiltonian, the energy is 
independent of Q.. We also suppose that, because the giant-resonance contribution to 
the deformation of a minimal-energy state arises from core polarisation of the nucleus 
by the quadrupole field of the SU(3) states, the minimal-energy state should be one for 
which the polarisation is maximally aligned with the SU(3) deformation. Inspection of 
Fig. 9 shows that such an alignment is achieved when the rotational angle Q! is zero. It 
then remains to determine the scale deformation for which E{d) is a minimum. 

We start from the following observations: 

(i) By construction, the lowest- weight state \o) = \No{Xoilo)) is an eigenstate of an 
isotropic harmonic-oscillator Hamiltonian 



Ha 



a2 I 



(48) 



with eigenvalue given by 

Ho\o)=Noh(o\o). (49) 

(ii) Every state \o{g)) E I^ is a lowest-weight state for an Sp(3,]R) irrep with respect 

to some harmonic-oscillator Hamiltonian. 
(iii) The lowest- weight state |(7(d)) is an eigenstate of a triaxial harmonic oscillator 



i?d = I 



2M 2 



^ + -M {wfxl, + coixi2 + «3<3) 



(50) 



We now determine the frequencies {cOt} of this triaxial oscillator for the minimal- 
energy state by use of the self-consistency property of mean-field theory which is 
that the minimum-energy state should be one for which (to the extent allowed by the 
mean-field constraints) the equidensity surfaces of the state have the same shape as the 
corresponding equipotential surfaces of the independent-particle Hamiltonian of which 
it is an eigenstate. 

It is important to understand that, although the lowest weight la) is observed to be an 
eigenstate of an independent-particle Hamiltonian, there is no implication that \o) is an 
independent-particle state. However, in many important cases, e.g., for an Sp(3, M) irrep 
built on a leading SU(3) irrep, it may be. 

An ellipsoidal surface associated with the density distribution of the lowest-weight 
state I o) is defined by the expectation values 

(x2)^ = (a|£x2,|a) = ^, /= 1,2,3, (51) 

and, unlike the isotropic harmonic-oscillator Hamiltonian (48) for which this state is an 
eigenstate, it is not spherically symmetric. On the other hand, a state \o{d)) which is an 
eigenstate of an anisotropic harmonic oscillator with potential 

Vix) = ^m£ {(ojxf,, + oixl2 + 0)143) . (52) 

has the mean values 

(^')a(d) = (a(d)|£4/l<7(d)) = ^, (53) 

An ellipsoidal equi-density surface for this state, defined by the equation 

2 2 2 

X^ X^ X'j 

+ T^ ^ T-?r — = '^onst. (54) 



(•^l)<7(d) (•^2)cT(d) (-^3)<T(d) 
is then given by 



Ni N2 N3 ^ ^ 



whereas an equi-potential surface for the potential (52) is given by 

cofxl + 0^X2 + 03^x\ = 1 . (56) 

Thus, for these surfaces to have the same shape, it is required that 

Ni (Oi = N2(i>2 = A^3C03 = k, (57) 

for some value of k. 

The volume of the ellipsoid (55) is 



3 \ (010)2 6)3/ 



(58) 



Thus, in order that this volume should be the same as that for coi = tt)2 = 0)3 = ft), we set 



and obtain 



ft)lft)30)3 = ft)^ (59) 



k^ =NiN2N^CD^. (60) 



This means that the scale parameters for the minimal-energy state should be such that 

{^l)a(d) _ (0 _ NiCO _ Ni 



{xj)a CO, k {NiN2N3y/^' 



(61) 



The above self-consistency method extends to all orders the methods used in Refs. 
[25, 26, 27] to estimate the strength x = hco/No to leading order in X/Nq and 11/ Nq, 
of the coupling constant for the effective quadrupole-quadrupole interaction in a model 
Hamiltonian of the form H = H^^ — ^XQ ' Q- 

Several important consequences follow from this shape-consistency result. 

(i) Coupling to giant resonance states (essentially) doubles the values of SU(3) 
quadrupole moments. 

This follows because quadrupole moments of the state \a) are given in the SU(3) 
model in a lowest-grade irrep, in terms of the oscillator unit a = y^Mco/h, by 

a^{{2xi-xl-xl))^ = 2Ni-N2-N3=2?^ + ^, (62) 

a^\/l{i4 -4))<y =N2-N3= Mo, (63) 

Therefore, for the shape-consistent state \o{d)), they are given by 

?, ? 9 ?s 2N}-N}-N} 

1— A/'S _ AJ-2 

^2/2/^2 ^2\ . _ ^^2 ^^3 (f^r^ 



An expansion in powers of Xq/Nq and /io/A^o then gives 

a2(2x?-xl-4)^(d)=2(2Ao + Mo) + (2Ao2 + 2AoMo-iUo^)/A^o + ---,(66) 

a^^{4~4)a{d) = 2^0 - (2AoAio - ^i)/No + • ■ • , (67) 

which, to leading order in Ao/A'o and Ho/Nq, imply that the effective charge of 
quadrupole moments and E2 transitions in the SU(3) model is predicted to be 

ggff = (2 + higher order terms .) x the bare nucleon charge, (68) 

where the bare nucleon charge for an SU(3) irrep would normally be Z/A times the 
proton charge, 
(ii) Coupling to the giant resonance excitations lowers the energies of deformed states. 

To estimate the magnitude of this lowering, one should compare the expectation 
values of {o\H\a) and {o{A)\H\o{A)) for a realistic nuclear Hamiltonian. How- 
ever, in the spirit of the shell model in which a first energy-ordering of states is 
given by their harmonic oscillator energies, we consider the energy of the state 
\o{d)) for the shape-consistent harmonic oscillator Hamiltonian H^i of Eq. (50). 
This gives 



■^N^,N2,Ni 



(a(d)|^d|<7(d)) =3M = 3(A^iA^2A^3)i/3f7a), (69) 



which is less than the energy of the lowest- weight state for the same irrep relative 
to the spherical harmonic oscillator Hamiltonian 



^A'o 



: Nohco = {Ni+N2+N3)h(0. (70) 



A consequence of this result is: 
(iii) The lowest-energy shell-model states of well-deformed nuclei are unlikely to lie in 
the lowest-energy spherical-harmonic oscillator shells. 

This claim is based on the observation that an Sp(3, R) irrep is defined uniquely by 
the quantum numbers Ni, N2, and A^3 for which 

Nq = Ni+N2+N3, Ao = iVi-A^2, ^o=N2=N3. (71) 

Moreover, all the mean-field states {\o{g)), ge Sp(3, R)} for a given Sp(3,R) irrep 
have common values for these quantum numbers and are of lowest- weight for the 
Sp(3, R) irrep with respect to some (generally triaxial) harmonic oscillator Hamil- 
tonian. We have argued that, with respect to any reasonable nuclear Hamiltonian, 
the lowest energy mean-field state should be one for which there is shape consis- 
tency between the density of the minimal-energy state and the harmonic oscilla- 
tor Hamiltonian for which the minimal-energy state is an eigenstate. A property 
of standard Hartree-Fock theory, on manifolds of Slater determinants, is that the 
minimal-energy state with respect to a given nuclear Hamiltonian is also of min- 
imal energy with respect to the locally defined Hartree-Fock independent-particle 



Hamiltonian. If we assume this property to hold for the above described Sp(3, R) 
mean-field theory, for which the analogue of the Hartree-Fock independent-particle 
Hamiltonian is the anisotropic harmonic oscillator Hamiltonian for which a given 
lowest-weight Sp(3, R) state is an eigenstate, the minimal energy state with respect 
to the given nuclear Hamiltonian should also be a minimal energy state with respect 
to the locally-defined harmonic-oscillator Hamiltonian. 

As noted above, the lowest-weight state |c7) of an Sp(3,R) irrep {No{Xoilo)) with 
respect to an isotropic harmonic oscillator Hamiltonian is an eigenstate of this 
Hamiltonian with energy Nohco. However, the shape-consistent lowest-weight state 
\o{d)) of the same volume is an eigenstate of an anisotropic harmonic oscillator 
of energy eigenvalue E!^^,No,N3 — ^{NiN2N3y'^ho), which is less thanNohco unless 
Ni,N2, and A^3 are all equal. Thus, it is appropriate to compare the minimal energies 
Eni,N2,Nj, for the possible Sp(3, R) irreps that are available to a given nucleus. 



6.3. Identifying the lowest-energy Sp(3, R) irreps 

The above shape-consistent mean-field results suggest that an energy ordering of 
the Sp(3,R) irreps available to a given nucleus is given by the values of Ej^/^^j^j^^j^/^. For 
the given nucleus, we can list the possible values of A^i, N2, and N3 together with the 
corresponding values of Ef^^^p^^^f^^. The first few irreps in order of increasing Epj^^pj^^N^ 
are shown for ^^C, ^^O and ^^^Er in Table 1. 



TABLE 1. Minimal mean-field energies EniN2N3 in units of hco for the lowest-energy irreps of 
'^C, '*'0, and '*''*Er as defined in the text. 
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The first remarkable result is that for the light nucleus ^^C the lowest three Sp(3, R) 
irreps are consistent with the states observed in this nucleus as described in an Sp(3, R) 
model framework by Dreyfuss et al. [28, 29] whose results are shown in Fig. 10. In 
particular, the highly deformed second excited state, known as the Hoyle state, which 
belongs to the (Ao,/io) = (12 0) irrep, falls well below its unperturbed energy of AhO) 
in the spherical harmonic oscillator shell model. A similar, result is found for ^^O, for 
which the lowest three irreps are consistent with those observed and described [27] in 
an Sp(3,R) model framework in Fig. 11. For ^^O, the first-excited state belongs to an 
irrep whose lowest- weight state would again be AfiOi above its unperturbed spherical 
harmonic oscillator and would lie outside of any normal shell-model valence space. 
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FIGURE 10. The low-energy spectrum of C energies calculated for the first three Sp(3, R) irreps 
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shown for C in Table 1 by Dreyfuss et al. [29]. (The figure is from their preprint with permission.) 
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FIGURE 11. The low-energy spectrum of ^^O energies calculated for a few Sp(3, R) irreps. (The figure 
isfromRef. [27].) 



For a range of Erbium and Tungsten isotopes, Jarrio et al. [24] estimated the most 
appropriate Sp(3, R) irreps needed to describe the ground-state rotational bands of these 
nuclei directly from experiments. Their method was simply to fit the E2 transition-rate 
data with an Sp(3, R) irrep in its rotor-model limit. They determined, for example, that 
the mean values of Aq and /io, for the low-energy states of ^^^Er are approximately 
(Afl) ~ 80 and (/io) ~ 12. They also determined the mean values of these quantum num- 
bers for Nilsson model states at the experimentally determined deformation and obtained 
the values (Aq) ~ 95 and (/io) ~ 12. Given that there is undoubtedly a mixing of Sp(3, R) 
irreps due to symmetry breaking interactions, their results are reasonably consistent with 



the dominant Sp(3, R) irrep contributing to the ground-state band of ^^^Er being the irrep 
with the next to lowest Ef^^i^^f^^ energy given for ^^^Er in Table 1. But they are inconsis- 
tent with the much less deformed irrep with A'o(Ao/io) = 814(30 8) and a slightly lower 
value of En^!\/^Nj, which is the most deformed irrep that could be obtained in the stan- 
dard spherical shell model. What is especially notable is that to obtain an Sp(3, R) irrep, 
with the large experimentally-observed deformation, it is necessary to go to unperturbed 
spherical-harmonic-oscillator shell-model states ofllhco above the lowest-energy irrep 
ofNo = SU. 



6.4. The wave functions of a symplectic model calculation 



Several algebraic symplectic model calculations have been made, where by algebraic 
we mean calculations with phenomenological Hamiltonians that are low-order polyno- 
mials in the Sp(3, R) Lie algebra. Such calculations, which explore the physical content 
of the symplectic model and its idiosynchracies in preparation for more detailed micro- 
scopic studies, became relatively straightforward following the development of vector 
coherent state methods [30, 31, 32] for calculating the matrix elements of the Sp(3, R) 
Lie algebra. 

Results of a calculation [23] of the energy levels of the ground-state rotational band 
of ^^^Er are shown in Fig. 8. The corresponding wave functions, shown in Fig. 12, are 
particularly revealing. They are very different from those calculated both with algebraic 
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FIGURE 12. Amplitudes for the wave functions of the 7 = 0, 6, and 10, Sp(3.R) model states of figure 
8 in a U(3) basis. The distribution of basis states includes states from 6 major harmonic-oscillator shells 
of positive parity that span a spherical harmonic-oscillator energy range of I2h(0. Note that the SU(3) 
amplitudes of all states are essentially independent of their angular momenta /. 



and microscopic interactions for light nuclei' which tend to have their dominant com- 
ponents in the lowest spherical harmonic-oscillator shells and contributions from higher 
shells that fall off rapidly over the space of a few shells. In contrast, the amplitude of 
the lowest-grade U(3) irrep, 826.5(78 0), labelled 0(78 0) in the figure, which is the 
left-most amplitude shown in the figure, is seen to represent a minuscule component 
of the states of the rotational band calculated. Given that the lowest-grade U(3) irrep 
826.5(78 0) is already at an excitation of l2ho) relative to a spherical harmonic oscil- 
lator shell-model Hamiltonian, the amplitudes shown in this figure are a clear indication 
of the irrelevance of shell-model calculations carried out in conventional shell-model 
spaces. However, from a consideration of Nils son model states for rotational bands in 
highly-deformed nuclei, this comes as no surprise. 



6.5. Some conclusions 

The above results lead to the following conclusions: 

(i) The partial ordering of Sp(3, R) irreps by the spherical harmonic oscillator ener- 
gies of their lowest-weight states, Nohco, is inappropriate for deformed nuclei. 

Recall that spin-orbit interactions also brings about a significant re-arrangment 
of the single-particle energy levels of the spherical harmonic oscillator and that 
pseudo-SU(3) [34] and pseudo-symplectic models [35] have been proposed to take 
this into account. The above considerations suggest that the re-ordering of the 
relevant shell-model configurations due to deformation is an even larger concern 
in strongly-deformed nuclei. Even so, the mixing of symplectic irreps by spin-orbit 
and other interactions, notably pairing interactions, has to be considered and is 
addressed briefly in the following section. Thus, the potentiality for combining the 
pseudo-symplectic techniques and the above methods is worth considering, 
(ii) The symplectic model in combination with self-consistent mean-field considera- 
tions demonstrates the importance of the quadrupole component of the mean field 
and leads to the expectation of coexisting states of different deformation already in 
the low-energy domain. 

The best known examples of this in light nuclei are the Hoyle state at 7.65 MeV in 
'^C and the first excited, 4-particle-4 hole, state in '^O, Many other examples of 
shape coexistence are given in the recent review of Heyde and Wood [22]. 
(iii) The dominance of highly deformed states in large open-shell nuclei, such as those 
of the rare-earth region, is an indication that, while the self-consistency formula, 
which orders Sp(3, R) irreps by E^^pj^i^^, provides a first estimate of the lowering of 
deformed states in heavy nuclei, it is probably an underestimate. 

It would appear to be an observed fact that highly-deformed irreps fall below 
the lowest-energy states of the spherical shell model in large doubly-open- shell 
nuclei. However, it unclear how far below they fall. This suggests the importance 



' See the several results reviewed in Ref. [33]. 



of seeking evidence of spherical excited states in otherwise deformed nuclei, which 
might correspond to closed sub-shell configurations (see [22]). It also suggests the 
importance of obtaining more reliable results from mean-field calculations in the 
spaces of Sp(3, R) irreps with meaningful nucleon-nucleon interactions. 

7. MICROSCOPIC CALCULATIONS IN SP(3,R) BASES 

The concern of this section is with the use of the subgroup chain 

Sp(3,R)DU(3)DSO(3), (72) 

to define shell-model coupling schemes for microscopic calculations. For a complete 
scheme, the basis states of this chain are combined with complementary spin-isospin 
states classified by the supermultiplet chain 

U(4)dSO(4)dSU(2)5xSU(2)^ (73) 

in an isospin formalism, or their neutron-proton counterparts. 

Enormous progress has been made by Dytrych et al. [36, 37] in the development of 
efficient shell-model calculations in Sp(3, R) D SU(3) bases and more general multi- 
shell SU(3) bases for no-core shell model calculations. These developments are of 
fundamental importance for showing that the many-nucleon theory of nuclear physics 
can be derived with nucleon interactions obtained from experiment, meson theory, and 
effective field theory (see lectures of Machleidt). 

It is to be hoped that many of the essential techniques developed in [36, 37] will 
also apply to heavy nuclei. The problem is the explosion of the dimensions of the no- 
core shell model approach with the number of spherical harmonic oscillator shells in 
the active space of the calculation, which Fig. 12, for example, shows can be be very 
large for heavy deformed nuclei. Moreover, as this overview indicates, they are expected 
to lie well above the lowest-energy harmonic-oscillator shells. Thus, we explore a 
complementary approach for applications to well-deformed states in rotational nuclei 
which reduces the dimensions of the Sp(3, R) spaces needed for converged results by 
several orders of magnitude. 

7.1. Calculations within a single Sp(3, R) irrep 

We consider a generator coordinate approach based on a procedure proposed by 
Filippov, Okhrimenko, Vasilevsky, Vassanji, and others [38, 39, 40, 41] and further 
developed by Carvalho, and Rowe [42, 43, 44]. 

This approach is founded on the observation that the real submanifold !^. of Sp(3, R) 
coherent states (44) spans the Hilbert space of the corresponding Sp(3, R) irrep. Thus, a 
basis for a truncated subspace is given by a discrete set of states 

lA^o(AoiUo)Ol) =^(^,)^(d,)^(^fe)|A^o(AoiUo)), (74) 



where |A'o(Ao/io)) is the lowest-weight state for the irrep. The notable fact, shown in 
Ref. [44], is that converged results are obtained using this method, to a given level of 
accuracy, with very few points when chosen in an optimal way. 

Good choices can be made following a computation of the mean field-energy 
E{Q.dQ.') defined by Eq. (47). For a rotationally invariant Hamiltonian, this energy is in- 
dependent of Q., i.e., E{Q.dQ/) = E{dQ.'). Taking linear combinations of states of many 
D.i then amounts to angular-momentum projection and the generation of states with 
good angular-momentum quantum numbers for which techniques have been developed, 
for example, by Cusson and Lee [45]. With only a single diagonal matrix d, namely that 
for which E{d) is a minimum and only the null rotation for Q.' E S0(3), it is found that 
rather good results are already obtained. However, the results, for a fixed value of d do 
not allow for centrifugal stretching effects corresponding to changes in the mean values 
of the elements of d with increasing angular momentum. Such effects are included by 
allowing d to vary with the angular momentum or by taking a small fixed set of {d/} 
about the mean-field minimum value d. (The preponderance of experimental evidence, 
based on yrast B(E2) values, suggests nearly constant d values.) Optimal {dj} sets can 
be determined along the lines considered by Carvalho et al. [44] in consideration of the 
energy functions E{A) and the range of values of d about the optimal value d for which 
(A/o(Ao/io)|^^(d)^(d)|A/o(Ao/io)) is non-negligible. Similar choices can be made for 

Applications of this approach to date obtain promising results with just a few {dy} 
and neglect of the Q' degree of freedom. Neglect of the Q' degree of freedom has 
an interesting physical significance. Note that an A^o(0 0) irrep, in which Ql does not 
contribute, contains the spectrum of U(3) irreps 

{{A^o + '^(«i-«2,"2-"3)}, n = 0,2,4,...}, (75) 

where ni,n2, and n-^ are even integers and n\>n2>n-i>0. Thus, as illustrated in Fig, 
5, the U(3) spectrum of states for a generic Sp(3,R) {Nq{Xq11q)) is given by the tensor 
products of the SU(3) irrep (Aq/Jo) with the U(3) irreps in the set (75). With suppression 
of the Ql degree of freedom, some of the irreps in this tensor product will be missing. 
However, no irrep of the stretched product form 

{NQ + n{'XQ + ni-n2,llQ + n2-m)}, n = 0,2,4,... (76) 

is omitted. But, which precisely which irreps are missing is not obvious to me at 
this time. This needs investigation and, if necessary, the Ql degree of freedom can be 
included. 

The technology for the efficient application of the above generator coordinate pro- 
cedure remains to be fully developed. So far, it has been applied to Sp(3,]R) irreps 
of maximal space symmetry; these are spaces of zero neutron and proton intrinsic 
spins. More generally, it is necessary to consider an irrep of the direct-product group 
Sp(3, R) X SU(4), in an isospin formalism, or the tensor product of two irreps of the 
direct product Sp(3, R) x SU(2), in a neutron-proton formalism. 



7.2. Mixed-representation calculations in an Sp(3, R) D U(3) basis and 

quasi-dynamical symmetry 

The development of efficient programs for the calculation of states of even-even nuclei 
within spaces of multiple Sp(3,R) irreps of maximum space symmetry (i.e., spaces of 
zero neutron and proton intrinsic spin) will undoubtedly be achieved in the near future. 
Exploratory calculations within such spaces are needed, starting with model interactions, 
to identify the nature of any problems that might arise, before the bigger task of elab- 
orating apposite techniques for mixed representations with non-zero intrinsic spins and 
microscopic interactions is tackled. This task would appear to be challenging but not 
insurmountable. Given that a relatively small number of optimally-chosen generator- 
coordinate basis states is needed, for each Sp(3,M) irrep, to obtain an acceptable level 
of converged solutions, and because rotational structure is much more clearly defined in 
heavy nuclei, it is possible that microscopic calculations in heavy deformed nuclei will 
eventually prove to be little, if any, more difficult than for the less deformed states of 
light nuclei. 

Calculations in spaces of mixed Sp(3, R) irreps are necessary to understand the influ- 
ences on rotational structures arising from various Sp(3, R) symmetry-breaking interac- 
tions, such as spin-orbit interactions and pairing interactions. Some years ago there was 
concern [46] that, while systems with a dynamical symmetry can be handled with rela- 
tive ease, systems with two or more competing dynamical symmetries would generally 
be intractable. However, it was also known that, if the dynamics associated with one of 
a competing pair of dynamical symmetries was slow and the other fast, there could be 
an adiabatic separation of the variables along the lines proposed by Bom and Oppen- 
heimer [47]. This concept was therefore incorporated into the language of group theory 
in terms of the precisely-defined mathematical concept of an embedded representation 
[48]. Such representations were subsequently observed to arise approximately in many 
situations of competing dynamical symmetries and were described as quasi-dynamical 
symmetries [49]. 

For a review of systems with competing dynamical symmetries see [50] . Such systems 
commonly exhibit a quantum phase transition between a phase in which one dynamical 
symmetry is dominant and another phase in which a competing dynamical symmetry is 
dominant. Moreover, the transition between the two phases typically becomes sharper 
for a many-particle system as the particle number is increased. It appears that such a 
dominant symmetry is almost invariably characteristic of a quasi-dynamical symmetry. 

The prototype of a quasi-dynamical symmetry is given by a system that has the 
dynamical symmetry of a rigid rotor with irreps defined by fixed intrinsic values of its 
quadrupole moments. When the rigid-rotor irreps are mixed by dynamical symmetry- 
breaking interactions, the system becomes a soft rotor with a distribution of intrinsic 
quadrupole moments. However, in the adiabatic limit, in which the centrifugal and 
Coriolis forces are weak, the intrinsic distribution remains constant over a range of low- 
energy angular-momentum states of its ground-state rotational band. In isolation, these 
states then have the same relative properties as those of a rigid-rotor; their non-rigid 
intrinsic structures are only revealed by generally small but non-zero transition matrix 
elements to states of excited bands. 



The following example [51] of an A'^ : 
as function of a control parameter a by 



48 particle system with a Hamiltonian, given 



H{a)=HQ + {l- a)Hsm + a^sus, 



(77) 



illustrates a close approach to a sharp phase transition from one phase to another as 
a is increased past a critical point. The low-energy spectrum of this Hamiltonian is 
shown, as a function of a in Fig. 13. When a = the Hamiltonian has an SU(2) C 
USp(6) dynamical symmetry and the system is in a superconducting vibrational phase. 
And, when a = 1 the Hamiltonian has an SU(3) C USp(6) dynamical symmetry and 
the system is in an adiabatic SU(3)-rotational phase (with low excitation energies). 
For intermediate values of a the two dynamical symmetries are in competition. The 
remarkable fact is that the spectrum is characteristic of the SU(2) phase for a < 0.55 
and characteristic of the SU(3) phase for a > 0.60. Moreover, as the particle number is 
increased the transition point becomes increasingly sharp. 
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FIGURE 13. The low-energy spectrum for the model Hamiltonian (77) as a function of a for a system 
of A^ = 48 particles. 

The corresponding wave functions for three values of a are shown in Fig. 14. What is 
notable is the sudden change in the character of the wave functions as a increases past 
its critical value of ~ 0.58 As a increases above the critical point, the SU(3) amplitudes 
show a reduced mixing of irreps but, more significantly, the mixed SU(3) structure of 
each state of the quasi-SU(3) rotational band rapidly becomes the same for each of the 
states in the band. It is as if all the SU(3) irreps were equivalent. It appears that only 
SU(3) irreps with similar quantum numbers are mixed strongly by symmetry-breaking 
interactions. 

The reduction in the mixing of basis states belonging to widely different irreps of a 
dominant dynamical symmetry is illustrated and can be understood in the ^^C and ^^O 
calculations mentioned above. For these nuclei, the mixing of the spherical harmonic 
oscillator lowest-energy configurations with those that form the highly deformed first- 
excited states is much suppressed by the fact that all matrix elements of one- and two- 
body operators between the low-lying Sp(3, R) irreps of the very different deformations. 
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FIGURE 14. Histograms of the SU(3) amplitudes for the lowest energy L = to 8 eigenstates of the 
Hamiltonian (77) for four values of a. 



shown in Table 1, are precisely zero. These irreps can only mix indirectly via higher- 
lying irreps. For example, the lowest-grade states of the excited irreps of ^^C and ^^O 
require the excitation of four particles relative to the ground-state irrep. Consequently, 
all matrix elements between states of the dominant irreps contributing to the ground and 
excited states are identically zero for a Hamiltonian with only two-body interactions. 

8. SUMMARY AND DISCUSSION 



I hope to have demonstrated in these talks that the well-developed group-theoretical and 
algebraic methods relating to symmetry provide important insights and powerful tools 
for investigating nuclear structure. In conjunction with geometry, analysis, and other 
standard tools, such as mean-field and coherent-state theory, they provide fundamental 
tools for studying the many interesting phenomena displayed in nuclear physics. 

A message I want to emphasise is the importance of determining the domain of valid- 
ity of a successful phenomenological model and of pursuing its relationship with other 
more sophisticated models, e.g., the many-nucleon theory of the nucleus. Investigating 
the limitations of a model is every bit as important as investigating its successes; it is 
probably more important. And establishing its microscopic foundations is what enables 
the model to contribute to the deeper understanding of the system it models. 

I have focused on the emergence of the symplectic model, exposed its roots in the 
Bohr collective model, and pursued its role in the development of a microscopic shell- 



model theory of rotational nuclei. However, I would not want to create the impression 
that the symplectic model stands alone. The primary pillars of nuclear structure remain. 
They include: the shell model, the collective model, the unified models, the pairing 
model, and the Hartree-Fock and Hartree-Bogolyubov mean-field theories. The sym- 
plectic model builds on these foundations and, in parallel with the algebraic pairing 
model, it defines a shell-model coupling scheme that enables one to make use of what 
one learns, from model fits to data, in the design of a more complete microscopic theory. 
Like the pairing model, which supplies the standard J J coupling scheme of the spherical 
shell model and provides a shell-model theory of singly-closed shell nuclei, the sym- 
plectic model supplies a coupling scheme that relates to shell-model bases expressed 
both in terms of spherical harmonic oscillator and anisotropic harmonic oscillator wave 
functions. Thus, as discussed in the text, it provides a shell-model theory of particular 
relevance to strongly deformed rotational states of nuclei. 

Some possibly novel perspectives have been presented. However, what is satisfying 
that the symplectic model embraces the views assembled from the vantage points of 
many developments and sets them in a new light. For example, an extension of the 
spherical shell-model with spherical harmonic oscillator wave functions to the Nilsson 
model [11] with anisotropic harmonic oscillator wave functions was made already in 
1955. However, the extension was designed only for the purpose of describing the states 
of added nucleons interacting with a deformed rotational core nucleus. Nevertheless, 
it was envisioned by many that the Nilsson model should somehow be extendable 
to a shell-model for deformed nuclei with non-spherical basis wave functions. The 
symplectic model, which relates many-particle Nilsson model states to lowest-weight 
states of Sp(3, R) irreps, makes this vision a reality. Moreover, as we have discussed, it 
leads to a practical generator coordinate approach to Sp(3, R) model calculations with 
microscopic interactions and approximations to mixed representation calculations with 
the inclusion of spin-orbit and other Sp(3, R) symmetry-breaking interactions. 

The symplectic model gives the deformed shell model a fundamental foundation. In 
particular, it resolves concerns about completeness and orthogonality. In the first place, 
states associated with different Sp(3,R) irreps and different angular momenta and spin 
are strictly orthogonal. Moreover, it is known that a geometric space of Sp(3, R) coherent 
states spans the Hilbert space for an Sp(3, R) irrep. Thus, the states of each such space 
are orthogonal to all other such spaces and, together, they span the whole many-nucleon 
Hilbert space. Thus, it is only necessary to orthogonalise basis states with common 
values for the associated quantum numbers. This is an enormous advantage over similar 
projected Hartree-Fock methods both computationally and conceptually; e.g., it retains a 
knowledge of what shell model spaces and Sp(3, R) irreps are included in a calculation. 

In conclusion, it should be emphasised that the construction of a shell-model theory 
of doubly open-shell nuclei is far from a minor topic in nuclear structure physics. 
Certainly doubly closed- and singly closed-shell nuclei have contributed a lot to the 
understanding of nuclear structure. However, such nuclei represent a small fraction of the 
nuclei available for study. Moreover, it is rapidly becoming apparent [22] that strongly 
deformed states, generally associated with doubly-open shell nuclei, are common among 
the excited states of singly and doubly closed-shell nuclei in which both the neutron and 
proton degrees become active. 
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